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A new approach to determine the value of the zero-temperature thermal exponent 9 in spin 
glasses is presented. It consists in describing the energy level spectrum in spin glasses only in 
terms of the properties of the lowest energy droplets and the lowest droplet exponents (LDEs) A; , (9; 
that describe the statistics of their sizes and gaps. We show how these LDEs yield the standard 
thermal exponent of droplet theory 9 through the relation, 9 — 9i + d\i. The present approach 
provides a new way to measure the thermal exponent 9 without any assumption about the correct 
procedure to generate typical low-lying excitations as is commonly done in many perturbation 
methods including domain wall calculations. To illustrate the usefulness of the method we present 
a detailed investigation of the properties of the lowest energy droplets in two-dimensional Gaussian 
Ising spin glasses. By independent measurements of both LDEs and an aspect-ratio analysis, we find 
6 (2d) ~ —0.46(1) < 9ow(2d) ~ —0.287 where 9ow is the thermal exponent obtained in domain-wall 
theory. We also discuss the origin of finite-volume corrections in the behavior of the LDE 9i and 
. relate them to the finite-volume corrections in the statistics of extreme values. Finally, we analyze 

£Sj ■ some geometrical properties of the lowest energy droplets finding results in agreement with those 

recently reported by Kawashima and Aoki 20 . All in all, we show that typical large-scale droplets are 
not probed by most of the present perturbation methods as they probably do not have a compact 
structure as has been recently suggested. We speculate that a multi-fractal scenario could be at the 
roots of the reported discrepancies on the value of the thermal exponent 9 in the two-dimensional 
. Gaussian Ising spin glass. 

o 

O 

Despite three decades of work in the field of spin glasses major issues related to their low-temperature behavior still 
remain unresolved 1 . Although important achievements have been obtained in the understanding of mean-field theory 2 
the appropriate treatment beyond mean-field to include short-range interactions is yet to be found. The absence of 
a successful analytical approach to deal with this problem corroborates the present state of our knowledge, often 
misguided by a non-accurate, if not confusing, interpretation of the numerical data. This situation has generated a 
hot debate about the correct physical interpretation of the available numerical data. Leaving aside the long-standing 
controversy whether replica symmetry breaking is or not a good description of the spin-glass phase 3 , there are still 
unresolved issues which are not as striking but evidence our ignorance about some fundamental questions. 

One among these problems is the correct value of the thermal exponent in two-dimensional (2d) Gaussian Ising spin 
glasses (GISG). This question has received attention from time to time during the last two decades but not enough to 
settle it definitively and explain the origin of some of the reported discrepancies. The study of the low-T properties of 
the 2d GISG starts with the work by McMillan who proposed 4 that thermal properties in spin glasses are determined 
by the scaling behavior of the typical largest excitations (commonly referred as droplets) present in the system. This 
idea has been further elaborated and extended to deal with equilibrium and dynamical properties of spin glasses in 
a scenario nowadays referred to as droplet model 6 . The low-T behavior in spin glasses is determined by a spectrum 
of large scale gapless droplets with typical length L and energy cost E ~ L s , 9 being the thermal exponent. As 
these droplets correspond to flipping some domains of spins (assumed to be compact clusters), the energy cost of 
these excitations arises from the set of unsatisfied bonds on their surface. The striking low-T behavior in spin glasses 
arises from multiple energy cancellations occurring at the surface of the droplet. These cancellations can be seen as 
the result of a competition between energy and entropy effects: as the droplet becomes progressively larger there are 
more available conformations for the surface to minimize the energy cost of the unsatisfied bonds. In the absence 
of cancellations one would expect 8 = (d — l)/2. However, as these cancellations are very important, the inequality 
< (d — l)/2 holds and 9 is by far less than the maximum value (d — l)/2. The value of the thermal exponent 
9 characterizes the low-T critical behavior as it is related to the correlation length exponent v where £ ~ T~ v by 
the identity v — —1/9. McMillan also used domain- wall renormalization group ideas to introduce a practical way 
to determine the leading energy cost of these low-lying large-scale excitations . The method consists in measuring 
the energy defect of a domain-wall spanning the whole system obtained by computing the change of the ground 
state energy when switching from periodic to anti periodic boundary conditions in one direction. Several works have 
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used McMillan's method to determine the value of 9 in two and three dimensions 7 ' 8 . Hereafter, in order to keep 
the discussion as clear as possible, we will denote by Odw the estimate of the exponent 9 obtained by domain-wall 
calculations. The initial value for 9dw reported by McMillan is Qdw = —0.281(5) for pretty modest lattice sizes 
L = 3 — 8. Recent numerical results with much more powerful algorithms have reached sizes L ~ 500 and confirmed 
the initial result with much larger accuracy 22,24,27 9p>w = —0.287(4). These studies would definitively close the 
problem if it were not by the existence of other alternative estimates of the exponent 9, largely consistent among 
them, which yield a quite different value 9 ~ —0.47(2). We will denote this estimate by 9tf as several of these 
methods use transfer matrix 18 . However, a word of caution is necessary here as the Monte Carlo method and other 
approaches that are not based on transfer matrix methods report values compatible with that estimate. For instance, 
Kawashima and Aoki used another method to estimate the stiffness exponent 20 . The idea is to generate a droplet 
inside a box of size L x L that includes a fixed central spin, with the following procedure. First, the ground state 
is found with a standard algorithms (we will denote it by the reference configuration). Afterwards, the spins at the 
boundaries of the box are fixed and the central spin is forced to flip respect to the reference configuration. The 
droplet of minimum energy that includes the central reversed spin and does not touch the boundaries is computed. 
The spanning length of the droplets generated in this way allows to define the fractal dimension of both the surface 
(or perimeter for the two dimensional case) and the volume. It is found that these minimum energy droplets have a 
fractal volume dimension smaller than 2 and the thermal exponent is 9 = —0.42(5) in agreement with results obtained 
from MC methods 26 and heuristic optimization algorithms 19 . A similar study of minimum energy clusters in the 
three-dimensional Edwards- Anderson model also reports evidence that 9p,w is an upper bound to the actual value of 
the thermal exponent 21 . 

The accuracy of previous estimates is poorer than the values obtained through the domain- wall method as they deal, 
in one way or another, with all possible excitations and not only with the calculation of ground state energies. More 
recently, another method has been used to estimate the value of 9. It consists in perturbing the original Hamiltonian 
Ho with a term eP, where V stands for the perturbation and e for its intensity. For example, P can be the overlap 
between the actual configuration and the ground state of the original Hamiltonian Ho- As e varies the new ground 
state of the total Hamiltonian H = Ho + eP remains unchanged until a certain value e — e c is reached where a excited 
energy level of Ho becomes the new ground state of H. The overlap between the old and the new ground states as 
well as the value of the shifting energy provoked by the perturbation links its energy cost E with its size providing 
another way to estimate 9. We will denote by 9p the estimate obtained in this way. This method has been recently 
used in the 2d GISG by Hartmann and Young 23 reporting the following value 9p « —0.31. Although slightly more 
negative than 6dw, both 9dw,9p appear to be statistically compatible. Yet more accurate estimates are needed to 
confirm whether 9p = 9p>w 

This last method and the domain-wall method have in common the same feature, i.e. they perturb the original 
Hamiltonian in one way or another to probe the characteristic energy of excitations that are supposed to be the 
typical ones that determine the low-T thermodynamic properties. In fact, the estimate 9dw can be considered as a 
particular example of 9p, where the perturbation consists in reversing all the bonds in one of the surfaces of the box. 
This raises the important question whether the different estimates of 9p, obtained by considering different class of 
perturbations, are different. The question is rather subtle as there are numerical indications that indeed this could 
be the case. For instance 8 , measurements of 9p where the perturbation is a uniform magnetic field yield a value 
9p = —0.48(1) compatible with the other competing set of values 9tf- 

How is that the value of the exponent 9p could depend on the type of perturbation? This is a very difficult question 
to answer as our present knowledge is inadequate. We can offer only speculative answers. Strong discrepancies 
among different types of perturbations could arise if a multi-fractal scenario governs the statistics of excitations in 
spin glasses. By definition, in all perturbation methods the probed large scale droplets are those which minimize the 
energy cost but constrained to maximize the value of the perturbation for the selected droplets. Therefore, among 
all possible large-scale low-lying droplets the perturbation method selectively probes those that maximally overlap 
with the perturbation. A dependence of the value 9p on a given class of perturbations could arise if the perturbation 
selectively probes one or another topological property of the droplet. This rather awkward multi-fractal scenario 
is not new in the field of disordered systems. Multifractality is known to be present in the localization problem in 
the strongly disordered regime. A multi-fractal scenario would imply the existence of different critical exponents at 
T = depending on the type of perturbation applied. On the other hand, the fact that the value estimated for 9p 
when the perturbation is a uniform magnetic field appears to be consistent with the value 9tf, suggests that maybe 
some types of perturbation can probe the relevant excitations while others may not. These good observables, which 
probe the typical excitations, could be called neutral observables in the same spirit as this term has been coined to 
describe observable dependences of the fluctuation-dissipation ratio (i.e. the effective temperature) in glassy systems. 
Concomitant, this "perturbation class dependence" issue is presently also debated in the different (but related to a 
certain degree) field of glassy dynamics. 

If the hypothetical multi-fractal scenario holds, then we must face the question about what is the correct procedure 
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to determine the thermal exponent 0. As determines the free energy cost of droplets, the natural answer is that 
is given by the lowest value among all possible estimates, 

6 = min r {O v } . (1) 

With the present available data this relation suggests that the estimate Otf is the correct value of the thermal 
exponent and that Odw as well as many other estimates Op are only upper bounds to the true value. 

The question we want to address in this paper is the following. Is it possible to devise a method that is alternative 
to current perturbation methods, in which excitations are not selectively probed by the perturbation, but selected 
only according to the correct balance between energy and entropy? The main purpose of this paper is to show that 
the analysis of the statistics of the first or lowest excitations gives a positive answer to this question. As we will 
see, the method we propose in this paper yields a consistent estimate of compatible with the value Otf, therefore 
supports the result that Odw and many other Op are only upper bounds to the actual value of 0. A preliminary 
account of these results has already appeared in 9 . 

The paper is divided as follows. Sec. II describes the basis of the lowest droplet approach and introduces the lowest 
droplet exponents. Sec. Ill shows the results obtained in the 2d GISG. Sec. IV analyzes a method to extract the 
value of the thermal exponent 0. Sec. V presents a more powerful method to extract the value of the lowest droplet 
exponents based on an aspect-ratio analysis. Sec. VI discusses the origin of the finite- volume corrections to the value 
of the lowest droplet exponent 0\ as a problem of corrections in the statistics of extreme values. Sec. VII analyzes some 
topological properties of the lowest droplets. Finally Sec. VIII presents the conclusions. There are also two technical 
appendixes: Appendix A presents the heuristic argument that Oi = —d for Gaussian spin glasses, and Appendix B 
explains the transfer matrix method we used to obtain the lowest droplets. 



II. BASIS OF THE LOWEST DROPLET APPROACH. 



The purpose of this work is to show an alternative approach to determine the low-T behavior of spin glasses by 
studying the size and energy spectrum of the lowest excitations by introducing two exponents (A; and Oi) needed to 
fully characterize the zero-temperature fixed point. All along the paper we will denote this exponents as lowest energy 
droplet exponents, or lowest droplet exponents in short, and that we will abbreviate as LDEs. The exponent A; is the 
most important one and describes the probability to find a large-scale lowest excitation spanning the whole system, 
while the exponent Oi describes the system-size dependence of the average energy cost of these lowest excitations. 

The underlying theoretical background of the approach is the following. To investigate the leading low-temperature 
behavior in spin glasses let us consider expectation values for moments of the order parameter by keeping only the 
ground state and the first or lowest excitation. This approach was introduced in 11 and can be shown to capture 
the low-temperature behavior at the leading order. The method that investigates the low-T properties based on a 
restricted analysis of the spectrum to the absolute lowest excitations has been also used for the study of the localized 
phase in the disordered Anderson model 12 . The present paper can be seen as the applicability of these ideas to the 
spin-glass case. At the end of the paper (see Sec. VIII) we will give reasons supporting the validity of our approach. 

To generate the spectrum of lowest excitations we consider the following procedure. Let us consider a set of M s 
samples and for each of them we determine both the configurations of the ground state and the lowest excitation. For 
a spin model the lowest excitation has v spins overturned with respect to the ground state (so the overlap between 
the ground and that excited state is q = 1 — 2v/V, V being the volume of the system) and with energy cost or gap 
E. It can be easily proved that the lowest excitation must be a connected cluster which we will generically call the 
lowest droplet. If v s and E(s) denote the volume and excitation energy of the lowest droplet for sample s, in the limit 
where Af 8 is sent to infinity, we can define the following joint probability distribution 

i ^ 

P(v,E) = - r ^S(v-v s )6(E-E(s)). (2) 

Using the Bayes theorem, this joint probability distribution can be written as P(v, E) = g v P v (E) where 

v_ 

2 /.OO 

= ! ; / dEP v (E) = l to. (3) 
v=l J ° 

g v is the probability to find a sample such that its lowest droplet has volume v and P v (E) is the conditioned probability 
for that droplet to have a gap equal to E. In what follows, we separately discuss the scaling behavior of both 
distributions g v ,P v (E). 
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Before continuing, and for sake of clarity, let us make an important digression about nomenclature. There are 
two volumes involved in the problem: the volume v of the lowest excitation and the volume V of the lattice. If not 
stated otherwise we will refer to the volume v as the size of the excitation while volume will generally refer to the 
lattice volume V. Thus, when we speak about finite-size excitations we usually refer to excitations with v finite, 
and finite-volume corrections (which we will sometimes abbreviate as FVC) will refer to the corrections affecting the 
distribution (2) due to the finite volume V of the lattice. 



A. The lowest droplet exponent A;. 

The simplest scenario for the size distribution of the lowest droplets is that all sizes occur with uniform probability. 
The normalization condition (3) imposes g v ~ 1/V. This situation is encountered in the Id GISG 6,11 with both free 
and periodic boundary conditions. However, in the most general situation, this does not hold and low energy droplets 
are found with a probability that depends on their size v. The simplest and most general way to incorporate such 
a dependence is to assume an ansatz solution for g v that factorizes into a power law A/V Xl+1 with A; > and a 
coefficient A = G(q) which depends only on the overlap q between the ground state and the lowest droplet, 

a (4) 

The behavior of G(q) can be guessed in both limits q — > 1 (the case g — > — 1 is equivalent in models with time-reversal 
symmetry which are those we are considering here) and q — > 0, 

G(q — > 0) — > constant (5) 
G(q ^ 1} -* (i-^TT • (6) 

The first relation describes the scaling behavior for the number of droplets whose size scales with the total volume 
of the system. As these can only depend on the volume V, G(0) must converge to a constant. The second relation 
is consequence of the fact that the number of droplets with finite size v cannot depend on V in the large V limit as 
these are not affected by the boundaries. On the other hand, the distribution of finite size droplets g v is self-similar 
as can be seen by inserting (6) in (4) and using the relation q = 1 — 2v/V. This yields g v ~ l/w Al+1 , the same relation 
as for the large scale limit (5) where gv ~ l/V Xl+1 . A simple expression that interpolates both limits is given by, 

G ^ = (^+(n^TT)- ( 7 ) 

Note however that, despite its simplicity, expression (7) is only an interpolation and the most we can say about G(q) 
concerns its asymptotic behaviors (5-6). 

The ansatz (4), applied only to large-scale excitations, was proposed in 11 . Note that although g v is defined for 
discrete volumes, in the limit V » 1, the values of q for consecutive droplet sizes v — > v + 1 become equally spaced 
by Aq = 2/V. Therefore, in the limit, l'>l, the function g(q) = ^g v becomes a continuous function if expressed in 
terms of the variable q instead of the integer variable v, 

9(1) = ^G(q) . (8) 

A word of caution is in order. Although (4) diverges for q = 1, leading apparently to a violation of the normalization 
condition (3) for g v , it must be emphasized that no excitation has q = 1 so there is a maximum cutoff value q* = 1— 2/V 
corresponding to one-spin excitations. For instance, if we insert (8) into the normalization condition for g(q) we get 
in the large V limit, 
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implying A; > as expected since otherwise the normalization would not be possible in the large-F limit. The 
divergent term (q — > 1) in (8) shows that for A; > one-spin excitations are the most large in number among the 
whole spectrum of sizes. In fact, from (4), g(l) ~ 0(1) 3> g(V/2) ~ 1/V X+1 , so the majority of excitations have a 
finite size. However the average excitation size 
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V 

v = Y, v ^^v^ 00 V 1 - x ' , (10) 

v = l 

diverges in the V — > oo limit and differs from the typical excitation volume Vt yp ~ 0(1). Relation (10) provides a way 
to measure the exponent A; alternative to the use of the scaling behavior (4) . 

B. The lowest droplet exponent 9i 

The analysis of the gap distribution P V (E) goes along the same lines as we did for the distribution g v , but with one 
important difference. As the gap E describes the lowest among all possible excitation energies, it has to scale in the 
same way for all droplet sizes independently on their size (and, in particular, whether these are finite-size or large scale 
droplets). This statement refers to a scenario which hereafter we will call "the random energy-size droplet scenario" 
(RESD scenario) to specifically indicate that the distribution of the lowest energies of droplets is independent of their 
size. Mathematically it can be expressed as, 

P V (E)=P(E) , V«. (11) 

In addition, we follow the standard droplet model and assume that the spectrum is gapless and defined by an exponent 
9i which describes the characteristic energy of the lowest droplets whatever their size or overlap q with the ground 
state. If the scaling function P v (E) is independent of v it follows immediately that the non-conditioned or size-averaged 
gap probability distribution 

P(E) = J2 9vPv(E) = J2 avP{E) = P{E) (12) 

V>1 V>1 

where we used (11) and the normalization condition (3) for g v . From now on, if not stated otherwise, we will always 
refer to the size-averaged probability distribution P(E) with the clear understanding that it coincides with any of the 
conditioned distributions P V (E). As the spectrum of lowest excitations is gapless, the normalized distribution P(E) 
has the following scaling behavior, 

p w = (13) 

We stress that the exponent 9i is completely different from the standard thermal exponent (see next section) as 
they describe totally different excitations. The thermal exponent 9 describes the energy-length relation for droplets 
typically excited at finite temperatures while the lowest energy exponent 9i describes the droplets that are separated 
by the smallest gap respectively to the ground state, so that, in general, 9i < 9. 

We will argue below in Sec. II C that 9i = —d for a generic class of spin-glass systems with coupling distributions 
with finite weight at zero gap. In addition, this relation will provide an alternative interpretation of the lower critical 
dimension in terms of the exponent A; introduced in Sec. II A describing the properties of the spectrum of sizes of the 
lowest droplets. 

C. The standard thermal exponent 6 

Now we want to show how the exponents A/ and 9i combine to give the usual scaling exponent 9 describing the energy 
cost of typical thermal excitations in droplet theory. There are several ways to show this result. For simplicity, here 
we exemplify this relation by analyzing the low-T behavior of the second moment of the spin-glass order parameter 
at the order linear in T by keeping only the first excitation. If <7{ CT ,r} — y S» a i T i denotes the overlap between two 
replicas (i.e. configurations of different systems with the same realization of quenched disorder), then the expectation 
value (q 2 ) can be written as follows 11 , 

o poo / rp \ 

<q i > = l-—Y,j o dEP(v,E)v(V -v)sech 2 ^— J, (14) 
where P(v, E) is given by (2). A low-temperature expansion of (14) 11 ' 9 up to linear order in T yields, 
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v 

y 

<q 2 > = I - —J2 9vPv(0MV - v) (15) 

which shows that the leading behavior is determined by both g v and the density of states at zero gap P v (0). In the 
standard droplet model, it is generally assumed that typical low energy droplets have an average size v = ^2 V vg v ~ V 
of the order of the system size (such as those generated by DW perturbation) and finite weight at zero gap Py (0) ~ 
1/L e where 9 is the thermal exponent. In principle, a single exponent 9 describes the scaling behavior of typical 
large-scale droplets with volume v oc V and determines the zero-temperature critical behavior. As these large-scale 
droplets are typical they occur with finite (therefore independent of V) probability gv ~ 0(1) while small scale 
droplets are simply irrelevant g v ~o(i) ~ 0- This yields, 

T 



<q 2 > = l-cjj, (16) 

where c is a non-universal stiffness constant related to the particular model. One of the most relevant results from 
the ansatz (4) is that both small and large scale excitations contribute to low-temperature properties. In general, let 
us consider any expression (such as (15)) involving a sum over all possible volume excitations. Restricting the sum to 
the large-scale droplets (v/V finite) the net contribution to such sum is proportional to VgvPv{0) oc L~ 6l ~ dXl V(0) 
(where V is the scaling function appearing in (13)). Coming back to (15) and using (4) and (13), we note that both 
small and large-scale excitations yield a contribution to (15) of the same order and given by, 



where c; is another constant (different from the constant c appearing in (16)). Identifying both relations (16) and 
(17) we obtain the general relation, 

e = 6i + d\i. (is) 

This relation shows how the value of 9 can be computed from A; and 9i. Through the study of a specific example, we 
will see later that the exponents 9\ and A; have strong finite-volume corrections arising from the corrections present 
in the statistics of the extreme values. However, we will present alternative routes to overcome this dependence and 
provide an accurate estimate of 9. 

Now we come back to the aforementioned argument at the end of Sec. II B claiming that in the large- volume limit 9i 
must converge to the value — d in the case of coupling distributions with finite gap at zero coupling. The details of the 
argument are shown in appendix A. The argument has two parts. First, it is proved that one-spin excitations provide 
an upper bound for the LDE 9i . Then it is argued that this upper bound holds also for any finite-size excitations (such 
as two-spin clusters, three-spin clusters, and so on). We will see below how this result is supported by the numerical 
analysis of the data. Let us also note that this result, in a RESD scenario (see Sec. II B) can be linked to the linear 
dependence of the specific at low temperatures, a result widely accepted, but that has been revisited recently in 10 to 
show that it has strong FVC due to the systematic FVC present in the value of 9\. Inserting 9[ = —d, (18) becomes, 

9 = rf(A ; - 1) . (19) 

This relation provides a way to distinguish the lower critical dimension did in terms of the average size distribution 
of the lowest droplets. According to (10) the relation \i(di c d) — 1 distinguishes a regime where the average size of the 
lowest droplet grows with the volume of the system to a regime where the average size of the lowest droplet is finite, 

d < d lcd : lim v(V) = oo ; A; < 1, 9 < (20) 

V^oo 

d > d lcd : lim v(V) = 0(1) ; A; > 1,9 > . (21) 

The marginal case A; = 1, 9 = is specially interesting as the average size v could be finite or diverge with the size 
but slower than a power law. This scenario corresponds to the mean-field behavior as replica-symmetry is broken 
in both the standard RSB 3 or in the TNT 17 (standing for trivial-non trivial) scenarios. Therefore, the study of the 
size spectrum of the lowest excitations in spin glasses can be very useful to find out the correct value of the thermal 
exponent in models without a finite- T transition (such as the 2d GISG) as well as establishing the correct low-T 
scenario in models with a finite-T transition. In the next section we apply all these ideas to evaluate the thermal 
exponent for the 2d GISG. 
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III. STATISTICS OF THE LOWEST ENERGY DROPLETS IN THE 2D GISG 



Several numerical works have recently searched for low-lying excitations in spin glasses using heuristic algorithms 13 . 
But, to our knowledge, no study has ever presented exact results about the statistics of lowest excitations. We have 
exactly computed ground states and lowest excitations in two-dimensional Gaussian spin glasses defined by 

H = — Jjj (Tj(Tj , (22) 

i<j 

where the are the spins (±1) and the Jij are quenched random variables extracted from a Gaussian distribution of 
zero mean and unit variance. These have been computed by using a transfer matrix method working in the spin basis. 
Representing each spins state by a weight and a graduation in the energy we can build explicitly the ground state by 
keeping the largest energy and, by subsequent iteration, the first excitation and so on (see Appendix B for the details 
on how we compute these quantities). The continuous values for the couplings assures that there is no accidental 
degeneracy in the system (apart from the trivial time-reversal symmetry a — > —a). Calculations have been done in 
systems with free boundary conditions in both directions (FF), periodic boundary conditions in both directions (PP) 
and free boundary conditions in one direction but periodic in the other (FP). In all cases we find the same qualitative 
and quantitative results indicating that we are seeing the correct critical behavior. 

We have found ground states and lowest droplets for systems ranging from L = 4 up to L = 11 for PP and 
up to L = 16 for FP and FF. The number of samples is very large, typically 10 6 for all sizes. The large number 
of samples assures us that many samples have large-scale droplets as first excitations. This provides us with good 
statistics to properly analyze the sector of large-scale excitations. The large number of samples requires a big amount 
of computational time so that calculations were done in a PC cluster during several months. For each sample we have 
evaluated the volume of the excitation v (and hence the overlap q = 1 — 2v/V between the ground state and the first 
excitation) and the gap E. From these quantities we can construct the g v and the P V (E). 




FIG. 1. g(q) versus 1 — q for the PP (left panel) and FF case (right panel) for different lattice sizes L = 5 — 11 (PP) and 
L=6-16 (FF) from top to bottom. In both insets we plot the scaling function g(q)V x versus 1 — q with A = 0.7. 



In Fig. 1 we show g(q) — ~g v as function of q for different sizes in the PP and FF cases. We can clearly see that there 
are excitations of all possible sizes but, as discussed in the paragraph following (9), the typical ones which dominate 
by far are single spin excitations. To have a rough idea of the number of rare samples giving large scale excitations 
let us say that nearly half of the total number of samples have one-spin lowest excitations, whereas less than 10% 
of the samples have lowest excitations with overlap q in the range — 0.5. This disparity increases systematically 
with size. For the lattice sizes explored the typical number of large-scale droplets is in the range 10 4 — 10 5 which is, 
indeed, quite good to have a good sampling of the sector corresponding to large scale excitations. A detailed analysis 
of the shape of g v reveals that it has a flat tail for large-scale excitations and a power-law divergence for finite-size 
excitations. The g v can be excellently fitted by the interpolating formula (7,8)), 
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(1-9) 



Ai + l 



(23) 



As shown in the insets of Fig. 1 a good collapse of the scaling function is obtained with the effective exponent A° ff ~ 0.7 
for both PP and FF cases. We also plot the line resulting from the fit of (23) with numerical data with the following 
values for A and B: PPBC: A = 1.55(3) and B = 0.777(3); FFBC: A = 2.02(3) and B = 0.85(1). Note that the fit is 
excellent and is hardly distinguishable from the points. The value of A; is compatible with the one obtained by fitting 
the average size with the expression (10) with the addition of a constant term to account for the small- behavior, 
v = C\ + CiV x ~ x . The same exponent A; can be estimated by measuring the ratio g(V/2)/g(l) ~ D\ + Z^V -1 ^. 
In both cases we get an effective exponent A^ ff = 0.70(5) as best fitting value. 

However these different estimates of A; arc strongly affected by finite- volume corrections (FVC). To evidence them 
we have estimated an effective L dependent X^ S (L) exponent by relating the average excitation size at consecutive 
sizes and using relation (10), 



Xf(L) = 1 - i 



(L) 



log(^ 



(24) 



In Fig. 2 we show A; (L) in the range L = 4 — 11 for the PP case. As we can appreciate there is a systematic increase 
of the effective exponent as we go to large volume sizes without any tendency to saturate. This proves that FVC in 
our measurements are still big and the estimate Af used to collapse the data in Fig. 1 is still far from the asymptotic 
exact value. 




FIG. 2. Effective lowest droplet exponent Af ff versus L for the PP case, computed using logarithmic derivatives. 



After having discussed the g v we jump now to discuss the scaling behavior of the energy gap distribution P V (E) 
and its average P(E). In Fig. 3 we show P(E) (main figure and inset a) and P V (E) (inset b) for the PP case. Similar 
results are obtained for the FF and FP cases. Quite remarkably, as was already anticipated in (11), the RESD scenario 
holds as the distribution P V (E) does not depend on the size v of the excitation (see inset b in Fig. 3), hence both 
large and finite-size excitations are described by the same gap distribution. 

In the main figure we can see how the width of distribution P(E) progressively shrinks to as L increases. Moreover, 
the P(E) has an exponential shape. This is shown in the inset a) of Fig. 3 where we plot P(E) in log-normal scale. 
Nonetheless, a detailed examination of the tails of P(E) reveals some deviations from linearity. In Sec. VI we discuss 
the origin of these deviations. We anticipate, though, that they are consequence of the strong FVC in the range of 
sizes investigated. In that inset we also verify the scaling ansatz (13) by showing the best data collapse for P(E) 
obtained with an effective exponent df ~ —1.7(1). This is very far from the expected value Q\ — —2 discussed in the 
preceding Sec. II C and in the Appendix A. A calculation of the moments of P(E) (13) for different values of L shows 
that there are also strong sub-dominant corrections to the leading scaling (13) that result in corrections as large as 
the ones affecting the exponent A;. 
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E 

FIG. 3. Gap distribution P{E) versus E for different lattice sizes in the PP case. In inset a) scaling obtained from the ansatz 
(13) with t9f ff = —1.7(1). In inset b) we show the P V (E) for different excitation sizes (q — 0.5, q = 0) for a lattice size L = 10. 
Note that the distribution is independent of the size of the excitation. 



Again, to manifest the magnitude of FVC in 9i we have evaluated E(L), the first moment of P(E), obtained by 
averaging the lowest gap over all possible droplet sizes for different lattice sizes in the range L = 4 — 11. We have 
estimated an effective L-dependent exponent by means of the following expression, 



9f(L) 



5 V E(L) J 



log 



L + l 
L 



(25) 



The results are shown in Fig. 4 for the PP case. Again, as for A® ff (see Fig. 2), we observe that the estimated value for 
6* z off systematically changes with size showing that, for the sizes we have explored, we are still far from the asymptotic 
regime. 
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FIG. 4. Effective droplet exponent 8f versus L for the PPBC case, computed using logarithmic derivatives (see text) 



We can summarize the results of this section saying that both lowest droplet exponents (LDEs) A/ and 0i display 
strong systematic finite- volume corrections (FVC). In principle, without further elaboration, it is difficult to give an 
accurate estimate for the thermal exponent 9 using (18). An alternative estimate for the exponent 9 could be defined 
from the analysis of the fraction of large-scale excitations with q < 1/2, f(q < 1/2) 14 , which is given by 







f(q < 1/2) ~ Vg(0) ~ l/V> 



L d{) 



l/L 6 



(26) 



where we have used 9i = —d (19). Although (26) yields estimates for 9, again these are affected by strong finite- 
volume corrections. In the range of sizes studied in this paper, and using (26) we get 9 ~ —0.6 quite far from the 
asymptotic value reported later in Sees. IV and V. How can we go further and estimate 9 in a safer way? In the next 
two sections we shall answer this question. 

IV. A GOOD ESTIMATE OF THE LOWEST DROPLET EXPONENTS 

An interesting aspect of the effective L-dependent exponents shown in Figs. 2 and 4 is that, while their FVC are 
large, their corrections are of opposite sign. While A° ff (L) increases with L, 9f s (L) decreases. As they have to be 
added to get 9 according to the relation (18) their finite-volume corrections cancel out to a certain degree. If we 
combine the two estimates for the best data collapse given in the previous section (A^ ff = 0.70(5), 9^ s ~ —1.7(1)) 
we obtain 9 ~ —0.3(2) which is very close to the DW value in average. However, this estimate is too pessimistic. A 
better route would be to use the two LDEs estimated from (24), (25) and adding them according to (18) 

9 cS (L) 



9f(L) + d\f(L) 



(27) 



In Fig. 5 (left panel) we show the value of 9 obtained in this way. Note that the value of the thermal exponent 9 has 
negligible FVC but relatively large statistical fluctuations with L. 

A better, albeit related, way to estimate 9 is the following. Instead of independently finding out A/ and 6>; we look 
for an estimator which depends on the appropriate combination of the two exponents 9 = 9i + dXi. The simplest 
quantity which satisfies this requirements is given by the combination, 



A(L) = L 



d E(L) 
v(L) 



(28) 



since E(L) ~ L 6 ' and v(L) ~ L d ^- X '\ using (18) we obtain A(L) ~ L e . To estimate the value of 9 we follow two 
different routes: 1) We use (25) by replacing 9f s (L) — > 9 cS (L) and E{L) — > A(L). By definition, this procedure gives 
exactly the estimate (27) shown in the left panel in Fig. 5. 2) A more stable estimate can be obtained from a fit of 
A(L) versus L, with data in the range [L, ■ ■ ■ , L max = 11] (for the PPBC case). This is shown in the left panel of 
Fig. 5 together with the previous estimate (27) and also in the right panel of Fig. 5 but there compared with the 
effective exponent 9dw obtained from domain-wall calculations. Our best value for 9 is 

9 = -0.46(1) . (29) 

This value is very close to the finite-temperature (Monte Carlo or transfer matrix) estimates 9tf = — 0.48(1) 19 but 
certainly smaller than the domain-wall value 9dw — — 0.285 7 ' 8 . Our estimate for 9 is compatible with the other- 
possible value 9tf obtained by other methods as discussed in Sec. I but is certainly inconsistent with the value 
obtained with other methods with results closer to the DW estimate. 
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FIG. 5. Exponent 9 for the PPBC case. Left plot: 9 exponent versus L obtained from two methods. Method 1: using (27). 
Method 2: using the more stable estimate fitting (28) over a given range of L values (see text). Right plot: Domain- wall 
exponent (top) and 9 exponent (bottom) estimated by the second method as explained in the text and plotted as a function of 
L. 
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All these estimates strongly support the inequality 9 — 9tf < &dw- However, one cannot exclude a situation where 
the present tendency of the data gets modified and 9 — » 9uw in the large- L limit 30 . We have already explained in 
Sec. IIC that 9i must converge to —2 in the large volume limit implying the relation (19). Introducing our estimate 
(29) in (19) we get, 

A; = 0.770(5) . (30) 

A convincing proof of the correctness of the values (29,30) requires proving that the estimate (24) converges to the 
value (30) when L —> oo. In the next section we present an aspect-ratio analysis to evidence that the estimates 
(29), (30) are correct in the large L limit. 



V. ASPECT-RATIO ANALYSIS OF THE LOWEST DROPLET EXPONENTS 

In this section, we present some additional data obtained via an aspect-ratio analysis (ARA). This analysis has 
been proved to be very useful to extract the value of the domain- wall exponent 9dw by generating domain walls in 
rectangular lattices M x L with different aspect ratios M/L 24,27 . It has been found that, in the limit of large aspect 
ratio, the value of 9ow for Gaussian spin glasses is largely independent of the boundary conditions. We have seen 
in Sec. Ill that our measurements on squared lattices of size L x L mix small excitations with large ones so one 
does not have a clear-cut separation in the statistical distribution between the two different regimes v ~ 0(1) and 
v/V ~ 0(1). Our main motivation here is to show that, by investigating large aspect-ratios, we can separate these 
two different scaling regimes. We made our measurements on systems of size L x M, with M — LR >> L where R 
ranges from 1 up to 10. We have investigated different types of boundary conditions: periodic boundary conditions 
in both directions (PPBC) and periodic boundary conditions in the L direction with free boundary conditions in the 
M direction (FPBC). 

In Fig. (6), we display the data for g(q) versus 1 — q (8) for the FPBC case for L = 8 and R = 1, 5 and 10. One can 
clearly see that the behavior of the distribution g(q) drastically changes as one increases R. Indeed, as we have already 
seen in Sec. Ill and in (5), (6), (8), for R = 1 it is very difficult to separate the region of small excitations (a scaling 
region with g(q) ~ j- — \ x +1 ) from the one of large excitations (a constant g-independent contribution g(q) ~ tt>t)- 

The main advantage of separating these two regions is that one can fit directly each of them. This yields two separate 
measurements of the LDE A; in addition to the estimate (10) obtained from the L dependence of the average size of 
the excitations. 
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In Fig. (7), we show g(q) versus 1 — q for R = 10 for various linear sizes L and for the PPBC case. These distributions 
have been obtained by running a large number of samples ranging from 10 million of samples for L = 4 down to 5 
million for the largest size L = 9. We have also inserted in the figure two vertical lines which indicate the limits for 
the range of values we have chosen for the fits of the scaling behavior of the finite-size excitation sector (l — q< 0.07) 
and for the constant contribution corresponding to large scale excitations (l—q> 0.25). We have chosen these values 
for the following reasons. First, as one can clearly see in Figs. (6,7), the scaling region for small excitations survives 
up to excitation sizes v ~ L x L. This size provides a threshold value for the overlap gth below which the simple 
scaling g(q) ~ ^_ g ^,+i does not hold anymore, 

2w N 2L 2 2 

1 - 9th = 1 - (1 ) ^ n ^ — ■ (31) 

yth v y) RL2 R \ ) 

An second, there is a crossover region around q ~ q t h- A careful look at Fig. (7) shows that the scaling region for 
small excitations ends around 1 — q ~ 0.07. At this value, one observes a change of the slope of the curves just before 
entering the regime of large excitations where g(q) becomes q independent. For 1 — q > 0.25, the curves are rather 
constant and the result of a fit does not depend much on the choice 1 — q = 0.25. This second threshold value is 
indicated as the rightmost vertical bar in Fig. (7). 
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FIG. 7. g(q) versus 1 - q for the PPBC case for R = 10. 

In Fig. (8), we show the estimated values of effective lowest droplet exponent A® obtained in three different ways. 
The first estimate has been obtained by averaging the volume of all excitations for different lattice sizes as explained 
in Sec. Ill and then taking a logarithmic derivative, see (24). The second estimate has been obtained by considering 
the large excitation sector (1 — q > 0.25) and its L, R dependence: 

g(q) ~ {RL 2 )~ X > . (32) 

Averaging the excitation volume within this sector (1 — q > 0.25) and using again the corresponding logarithmic 
derivatives as in (24) yields the second estimate. The third estimate for A; is obtained from a direct fit of g(q) for 
small values of 1 — q : 

g(q) ~ (1 - q)- 1 -* . (33) 

This third method is in fact the most direct one since it can be done for each size L (while the other two estimates 
require a fit using data from two different lattice sizes L and V). The first conclusion that we learn from Fig. (8) is that 
the ARA produces a great improvement on the estimated values of the exponent A; . The most stable measurement is 
the third estimate obtained by fitting the small-size spectrum of the excitations. In that case, Af ff is nearly constant 
with a value that converges to 
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A, = 0.77(1) 

in excellent agreement with the result (30) of the previous section. 
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FIG. 8. Effective lowest droplet exponent A^ ff versus L for the PPBC case for R = 10. We represent the values of X° s 
obtained from fitting the distribution of g(q) for small excitations (solid line), for large excitations (short dashed line) as well 
as the value obtained by fitting the average size of excitations (dotted line). 

Moreover, one also observes in Fig. (8) that the two other estimated values for Xf s , obtained with the first and second 
methods, are strongly correlated. This shows that finite-volume corrections, which are expected to affect the value of 
the exponent obtained from the analysis of large-size excitations, does affect also the value of the exponent obtained 
by averaging over the whole spectrum. In addition, we also observe that the ARA for large R strongly decreases the 
magnitude of finite-volume corrections. While on a square geometry, the effective exponent A® ff obtained from the 
average size of excitations took values in the range 0.52 — 0.62 (see Fig. (2)), with the ARA, we obtain for the same 
exponent values in the range 0.64 — 0.72, which are much closer to the expected asymptotic value 0.77(1). 
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FIG. 9. Effective exponent 6f 
curves which converge to 9f{L -> 



obtained via a logarithmic derivative for the PPBC and the FPBC. We also plot best fit 
oo) = -1.96(6) for the PPBC and to 6f(L -> oo) = -2.12(11) for the FPBC. 
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The same conclusion holds for the lowest droplet exponent 9[. In Fig. (9), we show the effective exponent 9f s 
obtained by evaluating the logarithmic derivative as in (25). Note that finite- volume corrections are much smaller 
than with the squared lattices and as a result, the value of the effective exponent converges much faster to the expected 
value -2. Using a fit of the form Of{L) = 9f{oo) + f£, one gets 6> ; cff (oc) = -1.96(6) for the PPBC, the best fit 
being also represented in Fig. (9). In this figure, we also show the same exponent obtained for the FPBC, where the 
best fit yields the asymptotic value 9^ s (oo) = —2.12(11). In both cases the fitting value we obtain for the exponent 
is a ~ 1. Note that the asymptotic values for 9f are well compatible with our prediction of Sec. II C, 9i = —2 (see 
also the heuristic argument in Appendix A). 

VI. FINITE- VOLUME CORRECTIONS (FVC) AND THEIR RELATION TO THE STATISTICS OF 

EXTREME VALUES 

What is the origin of these strong finite-volume corrections? Intuitively it is not difficult to find an explanation 
for the strong systematic finite- volume corrections in the lowest droplet exponent 0;. As the word lowest indicates, 
these exponents describe the statistical distribution of droplet excitations which are at the tail of the energy gap 
distribution that includes all possible high energy levels. As the volume of the system increases there is more available 
space to find excitations with lower energy gap. This implies that there is more probability to find a lowest droplet 
with an energy smaller than a given threshold value E* . As this probability systematically increases with L for all 
samples, 9f s (L) must be a decreasing function of L. This is in agreement with what we have found. The behavior of 
X° s (L) is more difficult to establish. By the same token, although there is more available space for droplets we expect 
that different droplet sizes increase or decrease their relative probability in a non-trivial way making difficult to guess 
how the exponent X[(L) systematically changes with L. 
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h 

FIG. 10. Local-field distribution for different lattice sizes with FFBC boundary conditions. 

To understand the origin of finite-volume corrections in the value of 9i we have focused our attention on the 
behavior of the upper bound exponent 9\ describing the statistics of the lowest one-spin excitations as described 
in the Appendix A. The gap distribution corresponding to these excitations can be obtained from the local-field 
distribution evaluated at the ground state. We have numerically computed this distribution for different sizes, the 
results are shown in Fig. 10. As discussed in the appendix A, the local-field distribution has a finite weight at zero 
field and is a self-averaging quantity. As the local field distribution is self-averaging, the probability distribution for 
the lowest one-spin excitations corresponds to the extreme value statistics of the local-field distribution p\ (h) where 
h stands for the local field which we assume to be positive as the gap is given by its absolute value (the subindex 1 
is used to stress that this distribution describes energy gaps for one-spin excitations only). If P\{h) stands for the 
probability distribution of the smallest local-fields, then Pi(h) can be easily related to pi(h) by standard probability 
arguments (see for instance, 28 ). Although the argument is very general, here we apply it to one-spin excitations. For 
a given sample, the lowest value h is selected as the minimum value among all the possible V local fields hi at each 
lattice site. The probability Pi(h) is given by the following expression, 



14 



P 1 (h) = V Pl (h)(l- ptWdh'y- 1 =-^\J } P(h')dtij (35) 

which accounts for all possible ways the value h coincides with the minimum value obtained among all different V 
local fields distributed according to the pi(h). The last identity shows that Pi(h) is normalized. This probability can 
be explicitly worked out in the large V limit, 

Pi(h) = -^exp[-Vgi(h)] = V gi {h)exp[-V gi (h)] . (36) 
Up to second order in h the function gi(h) is given by, 

From (36) we immediately learn that the gap distribution is an exponential with a sub-leading Gaussian correction 
whose magnitude decreases as 1/V. Actually, plotting Pi(h)/V as function of the scaling variable x = hV one gets, 



Pi(h) 
V 



2 



(38) 



In the large V limit g[(x/V) — > pi(0) and the coefficient in front of the Gaussian correction goes asymptotically to 
zero, therefore the distribution P\(h) converges to an exponential as expected, 

Pi{h) = Vpi(0)e X p(-Vpi(0)h) (39) 

in agreement with the scaling relation (13). We can now understand the deviations from the pure exponential behavior 
discussed in Sec. Ill in the context of the inset a) shown in Fig. 3. They are simply consequence of the finite- volume 
corrections of the extreme values of the gap distribution for all energy levels (and not only one-spin excitations as we 
are discussing here). So one could imagine to compute for a given sample the first V energy gaps corresponding to 
the first V excitations. We also assume that the resulting distribution p a \\(E) is self-averaging in the large volume 
limit (as it is the local-field distribution p\(h)). The lowest energy distribution constructed by taking the minimum 
value of the gap E for each sample yields the extremes distribution Pa\i(E) defined in (13) (from now on we will drop 
the subindex 'all' in P a \\(E) as it coincides with the P{E) defined in (13). Also the subindex 'all' for the Pan(E) 
will be dropped). In fact, the parameters p(0),p'(0) characterizing this distribution can be obtained from the P(E)'s 
shown in Fig. 3. To evaluate them, the best way is to analyze the cumulative distribution P(E) = J"J° dE'P(E') 
which from (35) we can assume to be P(E) = exp[— V g a \\{E)\. Thus we can fit V{E) with an exponential with 
Gaussian corrections Aexp{— Bx — Cx 2 /2] whose fitting parameters are related to p(0) and p'(Q). The best fits yield 
the following values p(0) « 0.2 and p'(0) = 0.3. 




15 



Coming back to our original goal we discuss now the finite- volume corrections for the estimate fff , as shown in 
Fig. 4. From the distribution (36) describing the whole spectrum of excitations we can express the effective exponent 
(25) for L » 1 as, 



The computation of E(L) is quite straightforward as it is given by the simple relation, 

E(L)=J^ EP(E)dE = J^ exp(-Vff(£0) (41) 

where we have used (36) plus an integration by parts. The integral, up to second order in 1/V yields, 

1 i/(0) + ( P (0)) 2 / 1 x 
Vp(0) V 2 { P {Q)f \V 3 J ' K ' 

Inserting this result in (40) we finally get, 



This shows that 9f (L) approaches —d from below (as p'(0) is positive). On the other hand the magnitude of the 

finite- volume corrections can be pretty large if q^o))^ >> 1- F° r instance, if one takes the results obtained from 

the analysis of one-spin excitations one gets pi(0) ~ 0.069, p[(0) — .125 yielding ^^yi — 27. which is indeed large. 

Inserting these values in (43) we obtain an estimation for 9f s (L = 12) = —1.65 in good agreement with numerical 
results (see Fig. 11). 

If we insert the previous estimated values for the whole spectrum of excitations extracted from the P(i?)'s in Fig. 3, 
we obtain ~ 7.5. From (43) it follows that 9f{L) ~ -2(1 - 7.5/V), which for L = 11 yields 9f = -1.87. All 

in all, the magnitude of the effective exponent 9i is well compatible with the reported value Of used in the inset a) 
in Fig. 3 for the PP case. Note that the FVC corrections to 9f s obtained from the local-field distributions in the FF 
case are much larger than FVC corrections in the PP case in agreement with ARA results (see Fig. 9). From this 
analysis it becomes clear that to significantly reduce the magnitude of the finite-volume corrections in the value of 9i 
(let us say 9i ~ —1.95), we would need larger volumes beyond 20 x 20. 



VII. COMPACTNESS OF THE LOWEST ENERGY DROPLETS. 



One intriguing question about the droplet excitations concerns their topological properties. Kawashima and Aoki 20 
have argued that droplet excitations are not compact. Instead, their volume has a fractal structure as the number of 
lattice points included in the droplet scales with its spanning length (which is a measure of the length scale of the 
droplet) with an exponent smaller (around 1.80(2)) than the dimension of the system (2). 

To answer this question we have computed the surface, i.e. the perimeter P, of all lowest droplets. The relation 
between the average perimeter as function of the size v of the excitation depends on both the fractal dimension of the 
surface or perimeter d s and the volume d v of the lowest droplets. d s and d v can be defined in terms of the spanning 
length I of the droplet which can be defined in different ways. For example, one could use the gyration radius, the 
average distance between the sites contained in the cluster, or the maximum distance among the sites of the cluster. 
As the typical length scale of our lowest droplets is small, I ~ 10, we have not attempted to estimate it as this 
can strongly depend on the precise definition of the spanning length. Here, we restrict ourselves to investigate the 
perimeter-volume dependence. In terms of the spanning length I the surface fractal and volume fractal dimensions 
d s , d v of the droplets are defined as, 

I ~ P~t (44) 
I - v^t (45) 

which combined give, 
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d 3 

P ~ v d -» 



(46) 



In Fig. 12 we show P(v) as a function of v for different lattice sizes. As can be seen, FVC are important for large 
volumes. However, there is an enveloping curve that is independent of L for small volumes and spans a progressively 
increase range of volumes as L increases. This enveloping curve is excellently fitted (continuous curve) by the scaling 
relation (46) and yields an estimate, 

~ 0.632(2) (47) 

d v 

consistent with the results reported by Kawashima and Aoki 4*- = 0.61(1) obtained with a completely different 
method. 
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FIG. 12. Perimeter (P) of the droplet versus its volume (v). The solid line corresponds to the fit (46) with d s /d v — 0.632(2). 



VIII. CONCLUSIONS 



We have shown that a proper description of low-temperature properties in two-dimensional Gaussian spin glasses 
can be done in terms of the lowest droplet exponents (LDEs) A; and 9i describing the spectrum of lowest excitations. 
A; describes the spectrum of sizes of the lowest energy droplets, while 9i describes the typical energy cost of these 
lowest droplets whatever their size. Assuming that 9i = —d one concludes that the LDE A; fully characterizes the 
spin-glass phase. Although independent numerical estimates of 9\ and A; show strong finite-volume corrections, the 
thermal exponent 9 = 6i + dXi can be well estimated giving the results (29), (30) 

9 = -0.46(1) A* = 0.770(5). (48) 

unambiguously showing that 6 < Qdw — — 0.287(4) 27 . Our estimates (48) have been confirmed via an aspect-ratio 
analysis which provides estimates much less influenced by finite-volume corrections. Moreover, the result 6i = —2 
(that is believed to be correct for spin glasses with coupling distributions with finite weight at zero coupling, see the 
Appendix A) has been numerically confirmed by the aspect-ratio analysis. To sum up, McMillan's excitations are not 
the typical low-lying excitations and our approach offers a new and independent way to estimate the thermal exponent 
9 without the need to generate typical low-lying excitations by looking at the new ground state of the system after 
perturbing it. 

We think that discrepancies on the value of the thermal exponent 9 reported by comparing non-perturbative methods 
(such as finite-temperature transfer-matrix calculations and the present lowest droplet analysis) with perturbative 
methods such as domain-wall calculations (or perturbations induced by introducing a coupling term in the energy 
function that induces a large-scale excitation) are serious enough to be taken as a clear indication that our knowledge 
of the low-temperature properties of the 2d GISG is still inadequate. In this direction we want also to recall the issue 
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of multifractality and the possibility that different exponents could describe the zero-temperature critical point. Is 
this really possible? Well, to our knowledge no exact result precludes this possibility and, although purely speculative 
at the present stage, one should seriously think about it. Altogether, the present analysis suggests that the excitations 
in 2d GISG are very different from the compact droplets proposed in the context of the droplet model. If this were 
true, the implications of the 2d studies in larger dimensions could be important. There are many routes that can be 
followed to understand better what is going on and the origin of this discrepancy Certainly, with the outstanding 
accuracy of present algorithms to compute ground states in 2d, it would be very interesting to revisit again the 
analysis of the statistics of the large scale excitations generated by imposing a uniform magnetic field. "Old" results 
by Rieger et al. 8 give an estimate for 9 that is compatible with our estimate rather than to the domain-wall estimate. 
This would be an independent check of our values, but using a perturbation method with an appropriate neutral 
observable such as the global magnetization as has been explained in Sec. I before (1). 

The proposed method may appear venturesome as, to our present knowledge, there is no numerical study in the field 
along this line of research. However, as explained in Sec. I, recent studies on the disordered Anderson model 12 have 
revealed that the analysis of the lowest excitation provides a good description of the localized phase. More studies 
are certainly required to understand better the reliability of the present method to investigate the critical properties 
of spin glasses. One disadvantage of our approach is that a huge number of samples is needed to reasonably sample 
large-scale excitations. However, as we saw in Sec. V, the behavior of the g(q) for small 1 — q can be extracted with 
a modest number of samples. The advantage, as has been already stressed in Sec. I, is that we do not introduce any 
external perturbation to generate the excitations. 

Finally, we want to comment on the extension of this approach to other models. Of course, the immediate extension 
one could think of is the 2d ±J model. However, the analysis of this model appears quite troublesome. This model 
does not have a continuous gap distribution but a discrete one that introduces further complications. As the ground 
state is not unique one has to redefine the full analysis to properly define the spectrum of lowest excitations. The 
discreteness of variables could have some unexpected effects in the present approach as seems to happen also with 
domain-wall calculations 22 ' 27 . It is more natural to extend the research to other models such as 2d ISG with other 
continuous coupling distributions without gap (e.g. characterized by P(J) ~ \ J\ a for |J| — > 0), Migdal-Kadanoff spin 
glasses (where both the ground state and the first excitation could be feasibly found with an appropriate algorithm), 
Gaussian spin glasses beyond d = 2 (where unfortunately, algorithms are much less effective than in 2d as the finding of 
the ground state becomes a NP complete problem) and finally mean-field spin-glass models where the zero temperature 
exponents are known and maybe the spectrum of lowest excitations could be analytically tackled. Preliminary results 
in this case 29 confirm that the present analysis describes pretty well the data for rather small sizes. We are pretty 
confident that, in the near future, new results and evidences will finally resolve this interesting problem. 

APPENDIX A: HEURISTIC PROOF OF THE IDENTITY 9 L = -D 

In this appendix we show that 9i = —d. In what follows we do not attempt to present a rigorous proof but we 
content ourselves to present an heuristic argument. The argument has two parts: first we show that — d is an upper 
bound, next we show that the upper bound is the exact value. For the upper bound the argument is well known 
and goes as follows. Consider the ground state and all possible one-spin excitations. Because one-spin excitations are 
not necessarily the absolute lowest ones, the statistics of the lowest one-spin excitations must yield an upper bound 
9] for the value of 9i, 6i < 9\. The statistics of the lowest one-spin excitations is determined by the behavior of 
the ground-state local field distribution p(h) in the limit h — > 0. If p(h) is self-averaging and p(0) is finite (in the 
large- volume limit) then the statistics of the lowest excitations must be governed by the exponent 6] = — d. Although 
we do not know a precise mathematical proof of the statement that p(0) is finite, it looks quite intuitive 16 . In any 
short-range system with a frustrated ground state and a coupling distribution with finite density at zero coupling, we 
may expect a finite probability to find a cage containing a spin coupled to its neighbors by a set of weak bonds which 
produce a vanishing net local- field acting on that spin. This argument should generally hold for d > 2. Moreover, as 
its name indicates, the local-field distribution is a local observable. An argument a la Brout proves that it should be 
self-averaging as all possible local field values are realized across the whole lattice 1 . The next part of the argument 
consists in proving that an identical upper bound is valid by considering excitations with size strictly larger than 
1 but finite. The upper bound derived for the one-spin excitations must necessarily hold for finite-size excitations 
beyond one-spin excitations (for instance two-spins, three-spins and so on) as the gap corresponding to the finite-size 



x Our numerical results in the 2d GISG confirm this conclusion, see Sec. VI 
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excitations can always be written as a linear combination of a finite number of local fields with coefficients which 
depend on the ground state configuration, ft is easy to verify that the aforementioned properties of the local-field 
distribution p(h) imply that the new gap distribution has a finite weight at zero gap and is self-averaging. This 
argument however cannot be extended to large-scale (with v ~ V) excitations in a straightforward way because the 
distribution for the corresponding gap distribution corresponds to an infinite sum of terms in the V — > oo limit. 
However once we argue that 9i is an upper bound valid for all finite-size excitations it can be concluded that this 
upper bound must coincide with the exponent Oi describing the probability of the absolute lowest excitations. From (4) 
the fraction of large scale excitations v — > V is given by Vgv = ttt • In general A > so this fraction vanishes 2 in 
the infinite- volume limit and finite-size excitations determine the result Oi — — d as they dominate the spectrum of 
lowest excitations. Moreover, if large-scale excitations yield a different value for Oi this would imply that boundary 
conditions could affect the value of the thermal exponent. That would be quite unusual as this would mean that the 
exponents of the T = fixed point would depend on the boundary conditions. 

APPENDIX B: TRANSFER MATRIX ALGORITHMS 

In this appendix, we will briefly explain how we determine the ground state and the first excited state. We will work 
on a square lattice of size L x L. The energy associate to a configuration of spins S(i,j) with a fixed configuration of 
disorder J x (i,j) and J v (i,j) is 

E= J2 E J x (iJ)S(i,j)S(i + l,j)+ E J y (iJ)S(i,j)S(i,j + l) (Bl) 

i=l,L-lj=l,L i=l,Lj=l,L-l 

+Bi J x (L,j)S(L,j)S(l,j) + B 2 J v (i,L)S(i,L)S(i,l), 

j=l,L i=l,L 

where B\ and B 2 correspond to the choice of boundary conditions. Here we will consider three cases : Periodic-Periodic 
boundary conditions (PPBC) with B\ = B 2 = 1; Free-Periodic boundary conditions (FPBC) with B\ = 1;B 2 = 
(or equivalently B\ — 0; Bi — 1) and Free- Free boundary conditions (FFBC) B\ = B2 = 0. We will only consider the 
case with a Gaussian distribution of the bond disorder J x , J y . To determine the ground state and the first excited 
states, we proceed as follow: we start by associating a weight for each configurations of spins in the first row of the 
lattice S(1,1),S(1,2), - ■■ ,S(1,L): 

W(S(l,l),S(l,2),.-.,S(l,L)) = B 2 jy(l,L)S(l,L)S(l,l)+ ]T J»(1,*)S(1,*)S(1,*+ 1) . (B2) 

1=1,1,-1 

Next, we start iterating the transfer matrix using a sparse-matrix factorization 31 . The first iteration gives 

^(5(1, 2), • • • , 5(1, L), 5(2, 1)) - max S (i,i) [•/*(!, 1)^(2, 1) + W(S(1, 1), • • • , 5(1, L))] . (B3) 

Since we are also interested in the first excited state, we define the second largest weight: 

H- 2 (5(l, 2), • • • , 5(1, L), 5(2, 1)) - mm s( i,i) [J x (l, 1)5(1, 1)5(2, 1) + W(S(1, 1), • • • , 5(1, L))} . (B4) 

In the following, we will use the simplified notation 

W(i,j) = W(S(i,j), • • • , S(i, L),S(i + 1, 1), • • • , S(i + l,j - 1)) . (B5) 

Thus eqs.(B3,B4) become 

W\l, 2) = max s{hl) [J K (1, 1)5(1, 1)5(2, 1) + W(l, 1)] (B6) 

and 

W 2 (l, 2) = rmn s{1A) [J x (l, 1)5(1, 1)5(2, 1) + W(l, 1)] . (B7) 



2 This fractions is finite only in d — 1 where A = 0. But this case is trivial as the surface of large scale droplets in d — 1 only 
contains a finite number of broken bonds. 
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At the next iteration, the numbers of possible weight will again be multiplied by two but we will keep only the two 
largest ones defined as 

W\l, 3) = max s(h2) [J x (l, 2)5(1, 2)5(2, 2) + J* (2, 1)5(2, 1)5(2, 2) + 1^(1, 2)] (B8) 



W 2 (l, 3) = max(min S (i,2)[J x {l, 2)5(1, 2)5(2, 2) + J y (2, 1)5(2, 1)5(2, 2) + W\\, 2)], (B9) 
maar S (i,2) [J x (l, 2)5(1, 2)5(2, 2) + jv(2, 1)5(2, 1)5(2, 2) + 2)]) . 

The general iteration relations are 

^(i.j) = maxs(ij-i)[J x (i, j - j - + 1, j - 1) (BIO) 
+ < 7"(i + l,j - 2)S(i + 1, j - 2)S(i + 1, j - 1) + W\i, j - 1)] 



W 2 (i, j) = max(min S ( i j-i)[J x (i,j - j - + 1, j - 1) (Bll) 

+J»(i + 1, j - 2)S(i + 1 , j — 2)S(i + 1, j - 1) + ^(i, j - 1)], 
max S (i lj -i)[J x {i,j - l)S(i,j - l)S(i + l,j - 1) 

+ 1 , j — 2)S(i + 1, j - 2)S(i + 1, j - 1) + W\i,j - 1)]) . 

In addition, each time that we end the construction of a new row, we must add the boundary term 

W\i, 1) -> W 1 ^, 1) + B 2 J v (i, L)S{i, L)S(i, 1) ; W 2 (i, 1) -» TT 2 (z, 1) + B 2 J y {i, L)S{i, L)S(i, 1) . (B12) 

We still have to take in account the boundary condition corresponding to B\ . The two type of boundary conditions 
(Free or Periodic) have to be considered separately : 

• Free boundary condition, B\ = : In that case, we iterate up to the construction of the weights associated 
to the configurations of the spins S(L, 1), • • • , S(L, L). The energy of the ground state (Eq) is then simply the 
maximum among all the weights W (L, 1): 

E = max {s{L , 1) ,.., S (L,L)}[W 1 (L, 1)]. (B13) 

We call {5°(L, 1), • • • , 5°(L, L)} the configuration of spins on the last row for the ground state. The energy of 
the first excited state is the second largest weight among I / F 1 (L, 1) and W 2 (L, 1): 

E x = max[max {s{L ^y... MLtL)} ^ { so (L ^ ) ^... tS o(L.L)}[W 1 {L 1 l)] 1 m . (B14) 



• Periodic boundary case, B\ = 1 : We first choose one configuration of spins on the first row 
5 l (l, 1), 5*(1, 2), • • • , 5 l (l, L). The weight of this configuration is defined as in (B2). The weight of all the 
other configurations of spins on this first row are fixed to an arbitrary large negative number. Next we iterate 
the transfer matrix as described above, up to the construction of the weights W X {L, 1) and W 2 (L, 1). Finally, 
we iterate one additional row, with bonds J X (L, i) and J V (L + 1, i) — 0. Next, we store the two weights W 1 and 
W 2 associated to the initial spins configuration 5*(1, 1), 5 l (l, 2), • • • , 5 l (l, L). We denote these two weights by 

W\S\l,l),S i {l,2),.-.,S i {l,L))= (B15) 
W\S(L + 1,1)- 5 l (l, 1), S(L + 1,2)- S\l, 2), • • • , S(L + 1, L) = 5 4 (1, L)) 

W 2 (S*(1,1),S*(1,2),---,S*(1,L))= (B16) 
W 2 {S(L + 1,1) = 5*(1, 1), S(L + 1,2) = S\l, 2), • • • , S(L + 1, L) = 5 4 (1, L)) . 



The energy of the ground state is the maximum on all the W 1 : 

E = max {S i {1 , 1) ,..,s^i,L)}[W 1 (S i (l, 1), 5^(1, 2), • • • , 5*(1, L))] , (B17) 

and we denote by 5°(1, 1), • • • , 5°(1, L) the configuration of spins on the first row for the ground state. The 
energy of the first excited state is the second largest weight among W 1 (i, 1) and W 2 (L, 1): 
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E 1 = raa a ;[m«a; {s , (M)i ... iS , (1;i)¥{SO(M)i ... ;S » (lii)} [W 1 (S i (l, 1), S\l, 2), • • • , S\l, L))], 

maa ;{S i (lil)i ... i si (liL)} [W 2 (S i (l ) 1), S% 2), • • • , L))]] . 



(B18) 
(B19) 



The construction of the ground state and of the first excited state is much more costly in computing time 
for the periodic case, since we have to repeat 2 L times the iterations, one time for each configuration 

So far, we have only described how to compute the value of the energies associated to the ground state and the first 
excited state. Since we also want to determine the spins configurations for these two states, we have to store, at each 
iteration of the transfer matrix, the value of the spin on which one sums, as well as the value of the previous spins. 
Thus, for each 2 L weights W 1 ^, j — 1), we have to store the configuration 

C S{1 ,,-i),.,s ( m ,,-2) = (S(l, 1), • • • , 5(1, L), S(2, 1), • • • , S(i, j - 2)) . (B20) 

At the next iteration, we will build the weight W 1 ^, j) with the corresponding configuration 

Cs(i,j),-,S(i+i,j-i) = (Cs(i,j-i),-,s(i+i,j-2), S(i, j - 1)) (B21) 

with S(i,j — 1) the value of the spin which corresponds to the maximum in (BIO). From this construction, one have 
access to the spins configurations of the ground state and the first excited state. 

Finally, one should also add that this construction can be easily extended to the second excited state, etc. After 
(B9), one can easily define a third weight which would be associated to the second excited state, and so on. 

Acknowledgments. We are grateful to J.-P. Bouchaud, D. Huse, O. Martin, M. A. Moore and A. P. Young for 
useful comments. We have been supported by the Spanish Ministerio de Ciencia y Tecnologia, projects BFM2001-3525 
(F.R.) and grant AP98-36523875 (M.S.). M. P. and F. R. acknowledge support from the French-Spanish collaboration 
(Picasso program and Acciones Integradas HF 1998-0097). Funding from the European Science Foundation through 
the SPHINX program is also acknowledged. 



1 A. P. Young, Ed. Spin Glasses and Random Fields, World Scientific, Singapore 1998. 

2 M. Mezard, G. Parisi and M. A. Virasoro, Spin Glass Theory and Beyond, World Scientific, Singapore 1987. 

3 E. Marinari, G. Parisi, F. Ricci-Tersenghi, J. J. Ruiz-Lorenzo and F. Zuliani, J. Stat. Phys. 98, 973 (2000); C. M. Newman 
and D. S. Stein, Preprint cond-mat/0105282. 

4 W. L. McMillan, J. Phys. C 17, 3179 (1984); A. J. Bray and M. A. Moore, J. Phys. C 17, L463 (1984). 

5 W. L. McMillan, Phys. Rev. B 29, 4026 (1984); ibid 30, 476 (1984); ibid 31, 342 (1985) 

6 A. J. Bray and M.A. Moore in Heidelberg Colloquium on Glassy Dynamics and Optimization, ed. L. van Hemmen and I. 
Morgenstern ( Springer- Verlag, Heidelberg, 1986); D. S. Fisher and D. A. Huse, Phys. Rev. B 38, 373 (1988); ibid, 386 (1988). 

7 A. J. Bray and M. A. Moore, Phys. Rev. B 31, 631 (1985); M. Palassini, A. P. Young, Phys. Rev. B 60, R9919 (1999); A. 
K. Hartmann, Phys. Rev. E 59, 84 (1999). 

8 H. Rieger, L. Santen, U. Blassum, M. Diehl and M. Jiinger, J. Phys. A (Math. Gen.) 29, 3939 (1996); 

9 M. Picco, F. Ritort and M. Sales, Statistics of lowest excitations in two-dimensional Gaussian spin glasses, Preprint cond- 
mat/0106554. 

10 M. Picco, F. Ritort and M. Sales, Random energy levels and low-temperature expansions for spin glasses, Preprint cond- 
mat/0106592. 

11 F. Ritort and M. Sales, J. Phys. A (Math. Gen.) 33, 6505 (2000); ibid J. Phys. A (Math. Gen.) 34 L333 (2001). 

12 Y. Gefen, R. Berkovits, I. V. Lerner and B. L. Altshuler, Phys. Rev. B 65, 081106(R) (2002); R.Berkovits, Y.Gefen, I.V.Lerner 
and B. L. Altshuler, to be published. 

13 F. Krzakala and O. C. Martin, Phys. Rev. Lett. 85, 3013 (2000); J. Houdayer, F. Krzakala and O. C. Martin, Eur. Phys. J. 
B 18, 467 (2000); A.K. Hartmann, Phys. Rev. E 63, 016106 (2001). 

14 This estimate was suggested to us by A. P. Young. 

15 A. A. Middleton, Phys. Rev. B 63 060202(R) (2001). 

16 We are indebted to D. Huse who pointed us this argument as well as the validity of the general result 8i = — d. 

17 M. Palassini and A. P. Young, Phys. Rev. Lett. 85, 3017 (2000); F. Krzakala and O. Martin, Phys. Rev. Lett. 85, 3013 
(2000). 

18 N. Kawashima, H. Hatano and M. Suzuki, J. Phys. A (Math. Gen.) 25 4985 (1992). 

19 N. Kawashima and M. Suzuki, J. Phys. A (Math. Gen.) 25, 1055 (1992); N. Kawashima, J.Phys. Soc. Jpn. 69, 987 (2000). 

20 N. Kawashima and T. Aoki, J. Phys. Soc. Jpn. 69, (2000) Supplement A 169-177. 

21 J. Lamarcq, J.-P. Bouchaud, O. C. Martin and M. Mezard, Europhys. Lett. 58, 321 (2002) 

22 A. K. Hartmann and A. P. Young, Phys. Rev. B 64, 180404 (2001). 



21 



A. K. Hartmann and A. P. Young, Large-Scale, Low-Energy Excitations in the Two-Dimensional Ising Spin Glass, Preprint 
cond-mat/0205659. 

A. C. Carter, A. J. Bray, M. A. Moore, Phys. Rev. Lett. 88, 077201 (2002). 
Kinzel and K. Binder, Phys. Rev. Lett. 50 , 1509 (1983). 

S. Liang, Phys. Rev. Lett 69, 2145 (1992); M. Ney-Nifle and A. P. Young, J. Phys. A (Math. Gen.) 30, 5311 (1997). 

A. K. Hartmann, A. J. Bray, A.C. Carter, M.A. Moore and A. P. Young, The stiffness exponent of two-dimensional Ising 

spin glasses for non-periodic boundary conditions using aspect-ratio scaling, Preprint cond-mat/0208546. 

M. Mezard and J. P. Bouchaud, J. Phys. A (Math. Gen.) 30 7997 (1997). 

M. Sales and F. Ritort (work in progress). 

M. A. Moore, Corrections to scaling in the droplet picture of spin glasses, Preprint cond-mat/0203469. 

M. P. Nightingale, in Finite Size Scaling and Numerical Simulations of Statistical Systems, edited by V. Privman (World 

Scientific, Singapore, 1990). 



22 



